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ABSTRACT 


This paper develops a model for computing the proba- 
bility of kill for an air defense complex composed of 
antiaircraft automatic weapons, radar controlled guns, and 
missile batteries. Two dimensional terrain was used to 
evaluate the model. The probabilities were determined at 
major terrain points along the route of approach to the 
Vital area for altitudes of up to 3000 feet above terrain. 
The curves of probability of kill versus altitude were 
found to be dependent on terrain, alr defense tactics, ame 
Weapon system parameters. A survivability index is calcu- 
lated by combining the probabilities of kill with a pilot 
Visual navigational probability. The resulting curves of 
survivability index versus altitude were found to be non- 
linear requiring a nonlinear programing technique to solve 
for the altitude of optimal survivability index within 
aircraft flight path constraints. The nonlinear solution 


was not included in this work. 
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CHAPTER I 


During the first year of the United States' involvement 
with the war in South East Asia, air power was employed 
extensively against many targets in North Vietnam. 6 
commanders of attack aviation units were required to plan 
and execute strikes against the same targets over and over 
again. Initially air defense of tnese targets was sparse 
to nonexistent. As time passed, air defense weapons were 
supplied from communist bloc countries until the problem 
presented to tne attack force commander was what to do 
about higher and higher attrition rates per raid. Specific- | 
ally, strike planners were seeking a solution to the problem 
Sr tinding the route of Approach and altitude to Tly= inte 
the target which would minimize aircraft losses from air 
defense means. One method of finding an answer to the 
problem is to analyze the losses from previous raids. This 
method is costly in lives and equipment and may not have 
produced data describing results from many altitudes and 
routes into the target. Using previous experience alone, 
one may never discover the optimal route. A major factor! 
in the effectiveness of land air defenses is the terrain. 
Very often terrain provides a natural route of approach 
Which in itself minimizes the capability of the air defense 
efforts. Sometimes this fact is obvious and in other 
situations it is not. It would be desirable in such κ. ^ 


situation to provide the attack air force commanders with a 


method that determines the optimal approach altitude and 
route that minimizes the air defense effort, The method 
Would generate an actual flight profile for the pilot. As 
new data on the enemy's air defenses are ascertained, up- 
dated flight profiles would be provided in the time period 
between strikes. The method referred to here is a digital 
computer model which evaluates all possible routes of 
approach into a target. For each route, the model evaluates 
the terrain, the air defenses situated therein, and the pilot 
navigational problem for all feasible altitudes up to a 
maximum value, and then determines the altitude along the 
route which maximizes aircraft survivability. with a 
functioning computer program which computes such a flight 
profile, new information could be inserted into the program 
and new results obtained in a matter of hours. Such a tool 
would be of tremendous value to attack air force commanders. 
Tne development of a mathematical model and computer 
solution for generating a flight profile which maximizes 
strike aircraft survivability was accomplished in a thesis by 
Lieutenant Colonel W. S. Miller, Jr., USMC, and Major E. E. 
Brown, USMC, reference l. In this work a model for calcu- 
lating probabilities of detection of strike aircraft over 
two dimensional terrain was combined with calculated proba- 
bilities of detection of ground navigational targets to 
generate a functional relationship of cost versus altitude 
of the flight path above terrain. The term cost is defined 


| 
to be a probability index which is a linear combination of | 


i 
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probability of radar detection and the probability that the 
ground navigational target was not detected at each alti- 
tude above terrain from 100 feet to 3,000 feet. Using this 
function as an objective function and calculating aircraft 
flight path constraints for each terrain point, the problem 
of determining that altitude which minimized the cost was 
solved as a linear program. The probability of detection 
mentioned above was calculated by computing at each terrain 
point and altitude the ratio of terrain visable to the 
total terrain within the coordinate system. As altitude 
increased, less terrain was masked from the point in question. 
NNeDefore, tne probability of detection was found to 1806265-6 
in approximately a linear fashion. This model assumed that 
a radar within the air defense system could be placed at 
any and all points along the terrain under the aircraft 
flight path. In other words, radar locations were considered ' 
wo be uniformly distributed over the entire route of approach, | 
To be sure, this method of calculating the probability of 
detection produces an indication of the effect of altitude 
on terrain masking in evaluating air defense capability. 
This paper extends the work of Miller and Brown. Since 
the uniform distribution of radar sites is not realistic and 
the method of calculating probability of radar detection is 
at best only an indication of air defense capability, a 
mathematical model and computer program will be developed to 
calculate the probability of kill for three types of air 


defense weapons sited in the terrain. Specifically, the 


mu 


problem to be solved will be to emplace an air defense com- 
ИС БЕС! ОВ "О: terrain in order to calculate atweach 
terrain point and altitude above terrain from 100 feet to 
3,000 feet along the aircraft flight path the cumulative 
probability of kill for an aircraft flying straight and 
level up to each terrain point. The terrain points con- 
sidered are those within the maximum effective range of the 
longest ranged air defense unit to the defended area. The 
problem reduces to finding the length of course line exposed 
to radar detection and within the effective range of the 
Weapon systems in order to determine the number of rounds 
that may be fired from which the probability of kill is 
determined, The model evaluates the terrain masking on 
detection range and firing time. The problem is developed 
for three dimensional terrain. The computer model was 
programmed for three dimensional terrain. No digitalized 
three dimensional terrain was found for use in this model. 
The computer program was accordingly verified using only 

two dimensional terrain which was available. That portion 
of the program which evaluates the three dimensional terrain 


Was not verified. 
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CHAPTER II 
THE MATHEMATICAL MODEL 


A mathematical model will be designed to compute the 
probability of kill for an approaching aircraft penetrating 
an air defense complex composed of heavy antiaircraft gun 
batteries, light antiaircraft automatic weapon batteries, and 
surface to air missile batteries. The purpose of this mese! 
is Lo determine a relationship of altitude versus probability 
at each major terrain point through the air defense complex. 
The model is designed to be general. No specific weapon ` 
system parameters in existence are used in the NES. E. 
Values for the parameters are of the order of those found 
in actual systems. The general nature of the model was 
selected in order to permit its use with any system and to 
prevent the necessity of classification under security 


regulations. 


PROBLEM FORMULATION 

Since it is necessary to compute the probability of 
kill at each terrain point and at altitudes of every 100 
feet above terrain to 3,000 feet, a concept for what is 
meant by a kill probability at a large number of specific 
points in space along a route of approach must be defined. 
The probability of kill by a weapon system is not defined at 
a point. Probability is calculated from the number of 


rounds the weapon system is capable of firing while the 
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target is within range, The kill probability at any partic- 
ular terrain point and altitude, PKxi> is defined to be the 
cumulative probability of kill resulting from the fire of all 
batteries of each type that bear on the target from each 
battery's maximum range to the terrain point in question. 

The amount of fire that each battery may produce is further 
constrained by terrain masking at the lower altitudes as it 
may occur, For example, at some terrain point k and alti- 
tude i, a line is drawn back along the route of approach 

to the beginning of the coordinate system or to the point 

of contact with a terrain formation, whichever is the shorter. 
Figure 1 is a diagram illustrating a two dimensional view 

of the terrain showing terrain points with altitude points 
above the terrain and the course lines associated with the 
(k,i)th point, representing the kth terrain point and ith 
altitude. In the figure, point F represents the 4th terrain 
point and the lst altitude. Drawn from this point is a 

line back along the course line of aircraft flight as pre- 
viously described. In this case, the line intersects a 
terrain formation at point G. The model assumes that when 
terrain interrupts a continuous flight, the aircraft will fly 
over the obstruction and continue on at the prescribed 
altitude. The line from G to F will represent a portion of 
the aircraft flight path. If a battery were located at 
point L, a probability of kill will be calculated from the 
number of rounds fired if sufficient time permits the firing 


of any rounds from any battery that bears on a target flying 
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the course from G to F, Such a battery would be located at 
pointe Ls. The value of probability calculated here will be a 
number greater than or equal to zero for altitude 1 at 
terrain point 4, These numbers become part of some func- 
tional relationship existing at terrain point 4. Another 
example of discontinuous flight is at points H and I. The 
Uninterrupted case 1S seen at point A; C; D and E: Here 

no terrain obstruction enters into the situation and the 
factors effecting tnie probability calculated depend on tne 
length of course lying in the detection radar ۳۷٦ 
and the weapon systeumweffectlveness envelope. As am eranp ios 
consider a battery located at terrain point 5. If the 

model were evaluating the probability of kill contributed 
by this battery up to point (4,2), the number of rounds that 
could be fired as the aircraft traveled from point B to point 
(4,2) would be the determining factor. Bis the point where 
the flight path passes into the radar beam pattern while 
point A is where the flight path enters the weapon system 
effedtiveness envelope, һе portici mof the fR in 

both radar and battery effectiveness envelope only, enables 
firing to occur, Since the altitude of approach is so low, 
the aircraft cannot first be detected until point B.. At 
point (4, 3) the battery at terrain point 5 detects the 
targets at point M in advance of its entry into the battery 
effectiveness envelope. A normal fire mission can be con- 


ducted in this case with adequate warning time. 
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Coordinate System 

In order for calculations to be carried out by thejcom- 
puter a coordinate system müst be placed on the tenrai 
section of interest. No existing map coordinate system 
will suffice. The coordinate square must be small enough to 
provide an acceptable approximation of the terrain É&ombours. 
Figure ረ is a diagram from a top view of the terrain and 
battery situation. This view shows two dimensions, distance 
along the course line, and offset distance from the course 
line for terrain and battery positions. 

The grid system shown in Figure 2 shall be the size 
necessary to cover the section of terrain holding the air 
defense complex of interest. A route of approach to be 
analyzed is selected and the grid system is placed over this 
route оп an appropriate map with tne x axis parallel to the 
@eurse direetien and centered ower the line. "At the inter- 
section of each grid line, the altitude taken from the 
Gerrain contour lines is recorded alomg with the X and Y 
Coordinates, The altitude above sea level becomes the H 
coordinate. All distances in this model are in feet. The 
X, Y, and H values are then punched on IBM cards for placing 
On cOmputer tape as data input. The job of preparing map 
data for transfer to computer tape is tedious and time 
consuming. For an individual to attempt the task would be 
impractical. A team of trained operations plotters equiped 
with the proper equipment would be capable of accomplishing 


the task. Since no three dimensional data was available, 
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the computer model was solved using two dimensional data. 
The model was verified with this cata except for the line 

of sight calculations for batteries not on the target flight 
line, The two dimensional terrain was taken from reference 
1 Each major terrain point was selected as тт a peak 
T a Valley. A sample of the terrain is shown in Figure 3 
and 4, Figure 3 is the detailed profile of the terrain. 
Figure 4 is the simplified version obtained by drawing 
straight lines between major terrain points. For a sample 
of the terrain actually used in the model see Figure 5 
тоц 7, Lae coordinates of the major peaks and valleys 
in tne xX direction and altitude in the H direction are re- 
corded in Appendix 2. The offset distance Y is zero for all 
computations but is carried along in the model for possible 
future use. Coordinates of the battery positions used are 
found in Appendix 2. Referring again to Figure 1, it should 
be noted tnat the X coordinate has its lowest value at the 
entrance of tne air defense complex and increases to the 


maximum X coordinate value at the defended area, 


Eur Defense Complex 

Tne three general types of weapons systems are the 
antiaircraft automatic weapon battery, the medium to heavy 
antiaircraft artillery battery, and the surface to air 
missile battery. These three are used since they cover the 
spectrum of ground to air weapon systems to be encountered 


in the reasonable future, 


1.2. 
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FIGURE 3 


SIMPLIFIED TERRAIN PROFILE 


FIGURE h 
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The automatic weapons battery is defined to be a 10 gun 
Ио арта firing weapons, five mounts, two guns to a 
mount. The fire control is by optics with target course 
and speed estimated by the crew, A computer is used to 
generate predicted azimuth and elevation for positioning 
the gun. Early warning is by radar with target acquisition 
by optical means. The medium to heavy antiaircraft artillery 
battery is composed of 4 guns with early warning radar, gun 
laying radar, and computer generated gun positioning signals 
via appropriate servo systems. 

The surface to air missile system is a low to medium 
altitude system. It is capable of detecting targets within 
its early radar beam pattern which is approximately from 1° 
elevation to 50° for volume coverage. The radar is assumed 
to possess effective moving target indication at the low 
altitudes. The tracking radars are also assumed to have 
low altitude capability. Since only one aircraft is attack- 
ing the complex in this model it is not necessary to specify 
the type of guidance with respect to the number of targets 
the system may attack at one time. 

The batteries are placed in the 2 dimensional terrain 
on the course line at various terrain points. The effect of 
terrain and increasing altitude, and the cumulative effect 
of more than one battery firing on a target is obtained from 
the various battery locations. No battery is permitted to 
fire at outbound targets, Each battery may fire from its 


maximum range or detection range up to its minimum range or 
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to its maximum turning rate in azimuth to the terrain point 
in question, whichever occurs first. No battery is permitted 
to fire on targets at a lower altitude than its own since 


this is not feasibly done. 


Weapon Target Kinematics 

In order to determine the length of target course line ` 
exposed to the fire of any battery, it is necessary бо μμ. 
the parameters of weapon target dynamics. Figure 5 48 a 
three dimensional diagram of the problem to be solved. The 
line C D is the path the aircraft flies as it passes into 
the range of battery J as evaluated at the (k, i)th point. 
Line A B is the projection of target course on a horizontal 
plane at the same altitude as the jth battery. The AH, the 
difference in altitude between the jth battery and the ith 


altitude of approach, is given by 


۸,,,, )ا ٭ کا 
AH [А መ : к=... 30‏ 
[Auk i By) | j = ውት EM : (1)‏ 


mae difference in position of the kth terrain point and the 


jth battery along the course line is given by 


Ж l vce N 
ል =|ቾፌ1” ጅሄ.)| Mee e 
ነ | у= l .... N, 


' (2) 


25 


m nm 
TARGET 


COURSE 






(k,i)th TERRAIN POINT 
C E i T, NX J 4 H r 





FIGURZ 8 


26 


Likewise for the offset distance of the jth battery 


= / N 
i = 9» 30 
AY - ITE ۱ _ |, N 
J ነ J == b. m. | እ 
(3) 
where 
( ==) | e € A 30 
= H + ioo! o TE 
AHK ; RU ዒላ = l, a 
(4) 
N = number of terrain points where a probability 


versus altitude function will be calculated 

30 = number of altitude points above terrain for 
which probability calculations will be accomplished 

Ny = number of Weapon batteries employed in 


ware terrain. 
The other parameters are defined as follows; 


ἃς 


(k,i)th point 


slant range from the battery to the 


ЕЕЕ z horizontal component of R 


R. 


s 


maximum effective range of the weapon system 


Bon = horizontal component of R, 
В, = horizontal component of the range at which 


the lower portion of the radar beam, at the estimated half 
power point of the radar beam pattern, intersects the target 


course line 
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X - the distance J L along the course line from 
the maximum open fire range of the system to the (k, i)th 
point being evaluated. This parameter is the key distance 
which will determine tne probability to be developed later, 


X =< Che warning distance, H.J from max"mam 


DR 
detection range to maximum open fire range 
Ө = the estimated elevation of the half power 


point of the lower edge of the radar beam pattern, 


Calculating Equations of the Model 

In this section all calculating equations will be 
explained or derived as appropriate. 
Prevent Outbound Firing. In a high threat operating con- 
dition when enemy attack is probable or imminent, air 
defense commanders are more concerned with the inbound raid 
than outbound raids. If a target passes through the fires 


of a particular battery, the probability of an outbound 


“kill is less in most cases than an inbound raid. The danger 


of firing at outbound raids results from not being ready to 
engage the next inbound xalq. ۲۲۲۶۶۰670۶ пе بی‎ oli o 
firing past the position of any battery is designed into the 
model. To insure no outbound firing AX is calculated in 


equation (2) without the absolute value notation. 


Aki Эх; ነ ሓጁ<ወ (5) 


m.s i | AX »0 
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Since any terrain point with a greater X coordinate value 
lies in an outbound direction from the jth battery, AX >0 

in this case. Firing is permitted up to the minimum range of 
the Jth battery, therefore setting AX = 0 permits later 


calculations for Xp up to Е= rather than to Ay ae Any 


J 
terrain whose X coordinate is less than Bx j Will produce a 
AX>0O and AX will be set equal to ИХ: - Bx;/. 

Determine Slant Range. The slant range to the (k, i)th 


point is given by 


- 2 
К. = Vax +A АНЕ (6) 


In Figure 8, it can be seen that R, is the vector sum of 

AX, AY, and AH. Therefore equation (6) will hold. For any 
battery j whose maximum range Ry < R, will not contribute any 
Probability of kill: sffierefore no calculations for the jth 
battery will be carried out in this case. И 
Determine if Target Elevation is Negative. Since firing at 
elevations below 0° is not feasible for many reasons, this 24 
model excludes such cases. After solving for AH in (1), 

if AH<O, excludes the jth battery from a probability 
calculation, Figure 9 illustrates this case. 

Minimum Firing Range. In certain cases which arise in this 
model the (k, i)th point lies over the jth battery or is 
inside the dead zone of battery effectiveness. The model 
must determine when this situation occurs and permit firing 


up to the minimum firing range for the ith altitude of 


approach. Figure 10 is a diagram of the minimum range 
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situation. Hg is the minimum firing range for the system. 
Angle B is the maximum elevation angle which is the minimum 
of the maximum radar antenna elevation angle and maximum 


elevation angle of the gun or launcher, 
Rxr = Re COS a (70) 
R as previously delfined is given by 


XS 
Rus = (R *- AH? (8) 


trea Figure ο. ΙΕ R... < Rep then firing must terminate at 


XS 
a horizontal range of Ryp rather than Hyg. Therefore equa- 


R {Ακ + ልክ” ) Rys > Rye 
° ar: (9) 


Maximum Azimuth Tracking Rates. When the jth battery is at 


tion (6) becomes 


some distance AY from the course line, it is possible for 
the maximum tracking rate of the system to be exceeded before 
the (K, i)th point is reached. If this should occur the 


distance X. would be decreased. The amount of this decrease 


D 
is calculated by finding that point on the course line where 
maximum tracking occurs and determining if it lies between 
(к, i)th point and the crossover point which is where AX = O 
or before the (k, i)th point is reached. If the point of 


maximum tracking rate occurs before (k, i)th point on the 


course line then Xp is reduced by the distance between the 


91 
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point of maximum tracking rate and (k, i)th point. Figure 
11 is a diagram of the azimuth tracking rate problem on the 
horizontal piane, Point K is the point where firing must 
terminate projected on the horizontal plane at the same 
altitude as battery J. The tracking rate must first be 
determined at this point. Let AZ be the azimuth tracking 


rate. Then 


A Z = Мт. 
К хз 1 
where У አጥ = tangential velocity Or the alreratft at Боза” 
and ቪ = horizontal component of the range to point K and 


Was calculated in equation (8). From Figure 11 it is clear 


tnat Von is given by 
۷7ھ‎ Ve AY 
AT Rus 


* 


Substituting this result into the equation above yields 


a ο 
AZ + ж» (10) 


4 


Now substituting equation (8) into (10) results in 
Az = GAT (11) 
RS--AM 
Dr AZ, = maximum tracking rate for the system engaged, and 
NZ < AZ. no constraint due to tracking rate is active. If 
however AZ > AZ,» then, the οσους το ας ге active апа fining 
must terminate at tne point where the maximum rate occurred. 


To find this point solve equation (10) for Ryg when AZ is 


p 


replaced by А2, resulting in 


tome — ሙ- = 


е. AY 
2 А “қа . (12) 


Using the value of R,, determined in (12) a new value for 


XS 
Rg is found by solving equation (8) for Rg as follows: 


К, = Re ant 


Equation (13) has changed the (k, i)th point to reflect the 


(13) 


constraint due to the maximum azimuth tracking rate, 

Maximum Elevation Tracking Hate. The same problem occurs 

in elevation as in azimuth for the possibility of exceeding 
tracking rates at low altitudes. In some cases the elevation 
tracking rate is exceeded at a range further from the (K, i)th 
| point than minimum range. It is necessary to compute the 
elevation tracking rate and determine if it is an active 
constraint. Figure 12 is a three dimensional diagram 
illustrating the vector component of target velocity tangen- 
tial to the battery target line of sight. The tangential 


velocity in elevation is given by 


. V, VAn* c Ay? 
Vere н Re (14) 


where the 


sin A = 


VAH*+ AY” 
ና 


3 
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Since the velocity Vyma 18 also opposite the angle A in the 


Smaller triangle 


XTE = Ут sin А 


Using the relationship as before that EL = 
М дн+ ду? 
> 1 
К. (15) 


XT 
“Rs”? results in 


E L = 


If ÉL» BL., where EL, ls the maximum tracking rate in eleva- 
tion of the particular system, then it is necessary to find 
the range at which the maximum rate occurs. Solving (16) 
for Rgp gives 


- መ ۷۷ ን ተ 
ы? 


(16) 


The new value of Rs is determined from 


and the new AX is determined by 
ጌ 2 
X * VR."—-AH-AY 
А Rs (17) 


Ate this point the "condita ons® fer terminating the firing 
at or before the (j, i)th point have been found. It is now 
necessary to find the point of open firing and then evaluate 
the effect of terrain on the length of course line between 
eae Open fire point, cease fire point, and=the battery \№0- 


cation. 
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Open Fire Point. Referring to Figure 8, it is seen that the 
OC Tire point is ar Hz The distance Ap is from the cease 
fire point which may or may not be the (k, i)th point to Н. 
Depending upon the value of H, the distance Xon may have 
positive, negative, or zero values, Point F must be found 

as it is the point where the detection radar first intersects 
Еле target course line., ΠΕΙΤΕ; diagrams the situation., 
Lines 1 and 2 represent different altitudes of approach and 
the different points of detection in range. The horizontal 


range to the point of detection is found as follows; 


Ah 
Rh = tan ወ 


(18) 
This relationship will be a good approximation since the 
ranges where the beam pattern begins to curve are far beyond 
the maximum effective range of any weapon system. The hori- 
zontal component of the maximum range of the system is 
determined to ре 


Rum = VR AP AY™ 


(19) 


The distance X) £ Хов on the target course line must be 
above ali terrain points and must be visible over its entire 
kength from the battery position in question. Te effect 

of terrain on low altitude approaches on the distance 

X) Ё Xon may be such that its value is decreased to less 


than zero when no firing is possible. A metnod of evaluat- 


ing the terrain's effect on the problem will be developed. 
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Terrain Evaluation. The method to be used for determining 
if line of sight exists from the point of radar detection to 
the battery position is to compute the slope of the line 
from this point to the battery as well as the slope from 
each intervening terrain point to the point of radar detec- 
tion. If the slopes of the terrain points are greater than 
the slope to the point of radar detection then line of sight 
exists. If one or more terrain points extend above the line 
joining the battery and the point of radar detection, tne 
line of sight does not exist for the length of the course 
line beyond the masking terrain points. A new point of 
radar detection is then found, A new slope line is com- 
puted from the masking terrain point at the altitude of 
approach and the battery. As before, all remaining inter- 
vening terrain point slopes are examined to determine if any 
more are masking the target course line. If another point 
has a slope less tnan the battery radar detection point 
line, then this higher terrain point becomes the new radar 
detection point. Depending on the terrain, the distance 

X) £ Ka may be uneffected or reduced to zero. If Xyx is 
the coordinate of the radar detection point then the slope 


to the battery position is given by 


g _ Αμ. = Өң; 
Lee Oj -Xak (20) 


The slope of each terrain point to radar detection point is 
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0 701 
2۴ Sed T Xue 
(21) 
where 
Алі = ат се 1. (1.4 ነ) ኮርጅ ርበብ ከ 
S = altitude of the (m,l)uh регла 7٣۳ 
m = mth terrain point between the battery and 


the radar detection point 

l = coordinate index of the offset distance 
inem the course line 

Xn 1 — X coordinate of the intervening terrain 
points, 

As long as the engaged battery is located under the 

target course line, the determination of major terrain 
points that may mask the line of sight is a relatively easy 
matter and tne slope calculations are simple. Figure 14 
is a diagram of the calculation of slope. Points 1, 2, 3, 
and 4 represent terrain points which may be typically found 
between battery 2 and the radar detection point L. As shown 
in the diagram, terrain 3 and ^ have greater slopes than 
tne line n They do not mask the line of sight. Terrain 
point 2 has a smaller slope than SOM and does mask tne line 
of sieht. Therefore, point 9 is the new radar detection 


point at which a new comparison slope must be calculated. 


Let Z be the distance from 43 to Вх: on the horizontal. Then 
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2 ፌሬ.” 65.1) (Cx X2) 
H 24 ы Šu, 


(22) 


This equation is derived from similar triangles in Figure 14. 


The value of the X coordinate of point 2 is given by 


КА = б, - 2. | (23) 


4hen the engaged battery is located at some offset 
distance, Y from the target course line, a procedure must be 
found to determine if line of sight exists between the 
battery j and all the terrain that lies between the lines 
Joining the cease fire point, and radar detection point on 
the target course line. Figure 15 illustrates this problem. 
Metice the coordinate grid. At each intersection of a grid 
coordinate, an IBM card recorded the X, Y grid index along 
with the altitude at the intersection. Notice also that 
tne X coordinate is parallel to the target course line. 
As previously indicated, the X coordinate increases in the 
direction of target flight. Point J is the battery position. 
me order to determine if line Of signt exists from J το Хук 
the slope of all terrain located on the line J, Xy aed LIL. ὃς 


computed and compared to the slope of the line J, X In 


NK” 
order to find the vertical and horizontal differential for 


all points on line J, X a series of triangles will be 


NK’ 
solved for the hypotenuse such as a, b, c in Figure 15. The 


distance ac is the horizontal distance for the slone 
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calculation m ihe distance bē is, if added To Ay? БИЕ ο 
coordinate of the closest point to P. Using the altitude 
of point P which has been recorded as previously described 
in digitalizing the terrain, the slope ОТ the line Xe” 15 


calculated, The distance be is given by 


be = ab tan a (24) 
where 

ab = distance between Y grid lines 

a 
tan a = XJ NK 
24 

us that 

ke = ab (Bx; - Жек) a 

АУ 

and 


ac = \ (ab)? + (be)* 


” (26) 
For each successive triangle, ab is replaced by Nab where 
Mae (B ድመ ) 3٦٢ 
y J к,1 
ала 


Бү, = У coordinate of battery j 


T: 1 = Y Coordinate of tne target course line. 
3 


3 


Substituting as appropriate, 


(Oxi =¥ ነ = e Bo > 
Y. ۱ (ob) + - 


oe 
іх = Ag ( By, = E (8x; 7 ይ i (28) 








Solving (27) and (28) for each value of N, enables the slope 


to be found at each Y. interval and establishes line of sight 


along the particular line being examined. If any terre 
slope is less than the comparison slope, then line of sight 


σος ВОВСЕ ОБЕ particular Line from J To Xy 1η 


NK' 


this case Хы is increased by one X coordinate grid distance, 


X, and the slope calculation process is again repeated. 


k 
The value ой Xa which remains after each slope line is 
compared with terrain determines the actual radar detection 
point. In the three dimens@enal case it Bis poseable vo 
track a target up to a terrain formation masking a portion 
of the target course line between the cease fire and radar 
detection point, and either continue tmacking 07 БИЕ 
mask in memory track, enabling relock, or to lose lock and 
require reacquisition on the other side of the mask. Іп 
the relock situation, the previous radar detection point 

X stands, If reacquisition is necessary, then a new 


NK 


Valhe of XNK is computed from the point where the terrain 


ceased to mask the course line, The value that Хук takes on 


ዓዱ 


is determined by 


τ- ἄν 


ዛሩ (29) 

aere 

I. = time for aircraft to pass terrain mask 

x. = distance on course line maskimg line of sight 
το the battery 

Va = aircraft speed 
nm pe Τη» where T is the memory tracking period, then Kx 
remains unchanged. If Ty > Tas Xak receives a new Value as 
shown below, 

Xak = Xg + Хк (30) 
where 

Ay = distance between X grid lines 
and 

N, = number of grid lines spanned by the terrain 
mask 


Determination of Length of Course Line for Firing. with 


the computation of K, and the cease fire point after con- 


κ 
Sidering maximum effective range of the system, maximum 
radar detection range, and effects of possible terrain 
masking, the distance Вх} - Ku 
after the system delays have been accounted for. The system 


is available for firing 


delays assumed in this model are as follows: 
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(1) Reaction time, ca which is the average of time 
for the radar operator to recognize target pips on the plan 
position indicator of the early warning radar set in the 
battery operations center, identification of the target and 
decision to engage the target. 

(2) Transfer action time, tas which is the time for the 
fire control operators to receive an assignment, lock the 
tracking radar on the target account for computer settling 
time, and gun crew reaction time to fire the first round. 

(3) Time of flight of the first round of the system to 
reach the target. 

Average values for (1) and (2) are assumed but time of flight 
is calculated from system parameters, An equation for 


solution of time of flight for gun systems is given by 


. - 
' Kp Kor (31) 
where 
Rof = Pange FO орет rire 
Ya < projectile muzzle velocity 
E = projectile drag 
If the k is known, (31) is a useful formula. If a system 


p 
graph of time of flight versus slant range is available, a 


regression equation is also useful. This model uses a re- 
gression equation approximateing a typical system curve of 
slant range versus time of flight, The concept is to place 


the first round on the target at the maximum effective range 


16 


of the system. If sufficient time exists forscht ur ee 
to travel from the point of detection to open Ότο ዓን በዓ 


then the system is assumed to expend the time 


tn = ty £ ta Z t. 


prior to the targets arrival at the open fire range, If 
the detection range is less than maximum effective range, 
енш ле time δη ub reduse Cae time for firing. If om 


nuc Other hand Only 3 part of the time Um Ts expended prior 


wo open Tire range, then some proportion of the ба will 
Geduce tne firing time. 

D, = tr Vr (32) 
Иеге 

D. спе свое весе travels during Une 


Ë 
system delay time 


AS a result Of terrain evaluar ron, X as Was seen, 


nk’ 

may nave a new value. 52۳۶۲٥۴6 the Tante to the maximum 

radar detection point is given by 
Rp = V(x ikl +A” (33) 

° 

This relationship is illustrated in Figure 15. Since detec- 

won range may have deoreased towbe Mess than 8) the open 


fire range, as a result of terrain restrictions, Hymn ls now 


even by 


LEN MIN AR Ram) (34) 
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Actually, Ben is the horizontal range component, The slant 


range to open fire becomes 





R = Rum + AX? + Ay” - 


Ihe distance along the course line from the radar detection 


point to open fire is given by 


X q ስመ " Rum = AX? ዩኒ δώρα 
OR (369 


Now the distance on the course line involved in reaction time 


δη may be all within XB if Xon > Uns Or Dari Of δη may be in 


XDR DI Ka < Uno If ቾ 58 = 0 then ሁጥ will reduce Xp by the 


amount сту Хр 


for the last round fired since any round fired after the 


Will also be reduced by the time of flight 


point (k, i) is reached does not contribute to the proba- 


51142 ዕዮዕ. 11. ነ ብ! ck) Therefore, Xp is given by 


πο μ᾽ - ልጆ-ዴ-ጫጭና ; ጾዕዚ*” 
X, ү О ےھ‎ AX-X tr Vr ነ Oc Xe 7) 
| Rim A? - ል።-ዣሻ ) ኣዕኗኞ Kr 
where 


T] т DR 


Ха = (oa 十 Vin 
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Having derived the equation for Ans the length of course 


line available for firing, it only remains to determine the 
number of rounds that will be fired in the time the target 


passes over X The time for target passage is 








D° 
т = А 
V (38) 
ፕ в 
Then the number of rounds fired is given by 
+ -- Ше 
ο E (32) 
ı/r 
where mE = time to fire one round and the rate of fire 
i 18 
8 t ° 
በ 
If Хов < Kas then the entire delay time including the time 


of flight of the first round occurs during the early warning 


time period and (39) becomes 


ДВ ; > 
Ж | "^pa 和 一 


+ $ Kor © Xr 


т 
Probability Model. In Chapter I, it was stated that the 


Ms 


(40) 





model was being developed for three different air defense 
systems. In the development of the model, the general 

scheme was to permit all calculations to apply to each type 
of weapon system except for the parameters used in the com- 
puter model, Except for the regression equations used for ` 
determining time of flight, this generality has been success- 
ful. However, the models must be separated for probability 


calculations into a separate one applying to each weapon 
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system. In the case of the two types of gun systems, a 
common model is still possible with one or two exceptions. 
The radar controlled antiaircraft gun system will be delt 
with first. The range at which each round is fired is 


computed as follows; 


К: = ско +ах- t а АНАУ? 


| (41) 


ας, бё, Ne 





This model assumes that the area a gun fires into is defined 
by a circular normal probability curve. The square of the 
radius of the circle described thusly has been determined 


empirically in references 2, 3, and 4 is given by 


2 = (105) +; * 15. (ET TY + (lot ves 


wnere 

o, = angular error of the gun in milliradians. 
gd; is the result of all the errors introduced by the fire 
control, data transmission, orientation, and alignment of 
the battery. Equation (42) is based upon the following 


theoretical equation: 


Б еи ^" — 
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where 
S = guasi static errors consisting of 
(1) Boresighting and alignment 
(2) Range errors (radar developed), radar 
timing errors due to mechanical alignment 
(3) Servo lag errors 
(4) Computer error 
(a) static errors 
(Б) ballistic match errors 
(c) prediction errors 
(5) Battery emplacement errors 
B+ 


"m 


Tracking perturbation 


[ 


Flight roughness 


The constants and results of using equation (42) have been 
verified in actual test firings according to reference 2, 
Target Area. The average cross sectional area of the tar- 
get's vulnerable sections is denoted by A in square feet. 

The vulnerable cross sectional area changes with different 
positional attitudes. An average value is assumed in the 
following calculations, In addition a method for calculating 
the area as a function of target attitude was also imple- 


mented, This area is computed from 


X 
Au Ard 


\ 


Ay + & А, ፪ Αμ (43) 


JL 


where 


A = percentage of area presented for the ith 
1 
round fired 


N = the range of the ith round 


2 
i 
at the ith round fired 


the X component of the range to the target 


Ay = cross sectional area of the standard target 


2 


defined to be 400 ft^ in the X coordinate. 


Ау = same as above for the Y coordinate set at 
800 ft* 

Аш = same as above for the H coordinate set at 
150 ft“ 

7 = the Y component of the range to the target 


at the ith round fired 
Н, = the H coordinate of the range to the target 
at the ith round fired 


The standard target may then be compared with any other 
target for converting probabilities. 

The probability that one round from the gun hits the target 
is the ratio of the vulnerable cross sectional area of the 


aircraft to the area the weapon fires into as follows; 


БЕС ee. = 
\= га 6 NR 
| ATO Шис (44) 


The probability required is the probability that one round 


kills the target. If the probability of kill given a hit 
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is PLK/H] then the probability of kill, assuming indepen- 
dence, is PLK/H] PLH]. Using this fact equation (44) 


becomes 


es —— PEKIN 


| aja” 
| (45) 


ШЕ: Na is the number of rounds fired, assuming independence 
of ronds fired from the gun, the probability of kill is 


given by 


Ne 
D one à Sa am 


=l 
where 
μα = number of guns in the battery. With all 
guns firing the probability that the aircraft is killed, 
assuming independence of guns in the battery, is given by 


ሠ = 1- (-ሚነ'ጭ 
۴۱ 
Е (47) 


The second type of weapon system is the automatic weapon 
battery. This battery does not possess gun laying radar. 
It depends on the optical means for acquiring the target. 
Early warning radar is assumed available. All equations in 
the previous model are applicable except equation (42). 
The total error for a gun with optical sights is assumed 
circular normal with the square of the radius of the circle 


of error given by 


g = (WV, te)“ (48) 


53 


This equation like (42) has been empirically determined and 
verified according to reference 2. 

In the surface to air missile battery, the third type 
of weapon system, a simplified model is assumed, Specific 
models are not feasible in order for the computer model to 
be applicable to any type of missile system, Therefore, the 
single shot kill probability is p. Then the battery kill 
probability is given by 

Q. = 1- (1-9) Αη 

J (49 ) 
Ka was calculated previously. However, in the missile 
battery the number of rounds fired at any target is con- 
strained by the basic load of missiles., Further the high 
single shot kill probability of a missile round precludes 
the firing of more than 2 or 3 missiles per target in most 
cases. This model insures that the value of مرا‎ permitted 


is as follows: 


ይፊ Ν΄» 2 
Мұс | ዙራ = | (20) 
а ከሩ! 


Any particular battery firing policy may change the value 
of N. as the tactical situation may dictate, 


For each terrain point k and altitude i, is cal- 


Pig 
culated for all batteries. Then the overall kill probability 


for all batteries is given by 
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ҚЫ. (51) 


ነ i ነ 
= ід sei, 3o 


The assumption of independence is quite valid for equation 
(51) since all batteries act independently in actual firing. 
Equation (46) also requires the assumption of independence 
between each round fired. If the weapon is a radar con- 
trolled gun using equation (42) to model the area within 
Which all rounds are fired, the probability of a hit for 
each round 79 و‎ must be independent for each value of i for 
the assumption to hold in (46). When examining the factors 
which determined (46), it is seen that certain of the errors 
are the same from round to round or are corrolated errors. 
These errors tend to break down the assumption of inde- 
pendence, Most of the errors are random in nature from 
round to round and are the factors supporting the assumption 
of independence. It is not known how serious the correlated 
errors effect the results of equation (46). If independence 
were not assumed, (46) would be modified as follows: 


[κ = Чазх -ዲ[5)(!-ዩ! P. e.) (1-6,0 ams اہ‎ 


p 


The difference between the probabilities in (^46) and the 
conditional probabilities in (52) constitutes the error if 
if is significant. All models found in the literature use 
the independence assumption. Investigation of the proba- 
bilities in (52) to determine how bad the assumption of 


independence really is may prove interesting. 
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CHAPTER III 
IHE COMPUTER MODEL 


In order to perform the calculations necessary to 
produce a relationship of altitude versus probability at 
each terrain point within the range of any battery in the 
defensive complex, a computer program capable of performing 
the calculations is essential. If the resulting computer 
program is to have any lasting value, the program itself 
must be documented in order to make clear how the program 
was designed and how the problem was solved by the computer, 
An explanation of the Fortran variables is found in Appendix 
III with the Fortran program located in Appendix IV. 

The program was originally designed to produce as an 
output an array of probabilities and corresponding altitudes 
for each terrain point examined. It was expected that the 
results would be punched on an output deck for use as an 
input to the computer program in reference 1. The output 
deck from this program would then replace the probability 
calculations in the other program, This replacement nas 
not been possible for reasons to be explained in a later 


chapter. 


Computer Inputs 


The inputs to the computer are of two types, The 
terrain digitalization and the weapon system parameters. 


The three arrays X d ic ord te 
array ail co an EM are the coordinates 


2n 


of the major terrain points being examined, In the two 
dimensional problem this is all the data needed. If the 
three dimensional data is available, three more arrays 
would be needed to record the terrain contours through out 
the area of interest. The arrays HH for altitude, XX 


ΤῊΝ ГУК 
X coordinate, and YY the Y coordinate of each recorded 


I,K 
terrain point would ай have toæbe inputs to the prognan: 

As previously explained, the model was developed for 
any weapon system. If the model is to represent any specific 
Weapon system the input parameters necessary to define the 
system must be clearly understood. An explanation of the 
input parameters is contained later in this chapter. The 


parametérs are defined in the Fortran code notation to 


facilitate their incorporation into the computer progran. 


Description of Computer Solution. 
Figures 16, 16a, and 17 are a block diagram of the com- 


puter program in general terms. The following description 
will relate to the block diagram. 

Block Zl. All data is read into the computer with units in 
feet and seconds,,angles are in degrees. 

Block #2. Most parameters of angular measure are given in 
degrees, The program converts all degrees to radiams. 
Block #3. The program is solved once for each N, batteries 
at each of 30 altītudes and for each of N terrain points, 


This constitutes 30 N.N iterations which in the sample used 


i 
in the actual program is 30 x 11 x 53 - 18,490 iterations. 
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Block Diagram of Computer Program 
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Block Diagram Continued 
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Block Diagram Subraitine LOS 
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CALCULATE TIE 
OF FLIGHT FOR 
GU™ ur 
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Figure 7 
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Block #4. In order to permit the program to use the same 
code for solving the problem for all three types of weapon 
systems, the system parameters are converted to a common set 
which represents the values corresponding to the particular 
system being examined at any one time, 
Blocks #5, 6, 7, 8, 9. Antiaircraft fire has been con- 
strained by 

(1) Firing only at inbound targets 

(2) Not firing at targets out of range 

(3) Not firing at targets lower in elevation 
than the battery 

(4) Not firing at less than minimun range 

(5) Not exceeding maximum azimuth and elevation 
tracking rates. 
A section of computer code as shown is Figure 16 provides the 
calculations to determine, at each iteration, the occurance 
of any one of the five constraints listed above, The overall 


effect of each of the constraints will now be explained. 


(1) The program is designed to permit firing up to 
minimum range or to the point of maximum azimuth and ele- 
vation tracking rates. If & terrain point is well behind 
the battery or in the opposite end of tne course from the 
direction of target attack, the minimum point is thusly 
established. On the other hand, if the terrain point in 
question is in the direction of attack, then the terrain 
point itself is the cease fire point unless the above 


constraints act to increase the cease fire range. 
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(2) Firing at targets that will always be out of range 
ls prevented by assigning zero to the probability of kill if 
the range to tne terrain point and altitude is beyond system 
range. 

(3) If a target is lower in elevation than the battery, 
zero is assigned the probability of kill for that battery 
at that altitude, 

(^) 4hen either minimum range, maximum azimuth and 
elevation tracking rates are exceeded, whichever occurs 
Burst, the cease fire range is increased to that point., 

lock #10. The subroutine LOS is called at this point which 
calculates the distance of course line available for firing. 
The function of LOS will be explained later. 

Block #11. The length of target course line available for 
firing resulting from terrain considerations may produce a 
new value for open fire range. Based upon the value of Хр, 
the open fire range is calculated. 

Block #12. Depending on what type of battery is engaged, 
the program will go through one of the three paths in this 
block, The probability of kill depending on the number of 
rounds fired is calculated for one of the three systems. 
Block #14, The overall kill probability is computed here 
for one iteration which includes all weapons at one terrain 
point and one altitude. 

Block #15. After all iterations are completed the program 
is capable of plotting the curves of altitude versus 


probability. Any other output is available when desired. 
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Subroutine LOS 

Block #1. Receiving all necessary data from the main pro- 
gram, all calculations planned are possible. Radar warning 
time is determined prior to terrain evaluation, The value 
for the point of first detection is the result of this 
biock, 

Block Z2. Here the program determines which terrain points 
lie between the battery and the point of first detection. 

If any terrain interrupts the line of sight, a new value 

for radar detection point is calculated. The computations 
are for two dimensional terrain which mean terrain for 
batteries on the target line of sight. 

Block #3, 4 If three dimensional terrain is used, sub- 
routine SEE is called where the point of radar detection, as 
a result of terrain, is found. 

Block #5. Open fire range without regard to delay is 
calculated for use in determining time of flight. 

Block #6, Again the program branches into three paths, one 
for each weapon system, This is necessary since each system 
has a different function of slant range versus time of flight. 
Block F . The delay time which is the sum of crew reaction 
time, system delay time, and time of flight of the first 
round fired is calculated in this section. This value is 
not used in the case where early warning time is less than 
total delay time. 

Block #8. After all the above calculations are completed, 


it is possible to find the actual point of open fire. This 


64 


point is the range at which the first round is placed during 
delay time to maximize the number of rounds fired, Three 
possibie situations arise in calculations for the open fire 
point as shown in Chapter II, equation (35). These are 

(1) Adequate warning time exists for ali delay time 
to occur during the early warning period, 

(2) Only part of the delay time occurs during the early 
Warning period and the open fire point is reduced, 

(3) All of the delay time must fall during possible 
firing time due to the radar detection point being less 
than maximum effective range. The subroutine sends X 


D 
back to the main program. 


A Description of Input Parameters 


At some later time, it may be desirable to compute 
the results from the model for an actual air defense system. 
In order to do this, an explanation of the input parameters 
is essential. A description of how to insert new inputs 
is provided below. Statement number 4 in Appendix III is 
the data statement, Contained in this statement are all 
the system parameters, The program is convertable to any 
specific weapon systems by changing the values in the pro- 
gram data statement to conform to those of the new systems. 
Included with the description of each parameter is the value 


assumed in the current form of the progran. 
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Parameter 


NAP 


NEP 


NFU 


ELRMG 


ELRMA 


ELRMM 


AZRMG 


AZRMA 


AZRMM 


RMG 


THERI 
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Description 


Number of altitude points 
above terrain examined, 
ΝΑΡ - 30 


Number of ሞጋ points 
examined. NAP = 84 


Number of batteries placed 
in the terrain. NEU = 11 


Maximum elevation tracking 
rate for the radar controlled 
gun system. ELRMG = 10°/sec 


Maximum elevation tracking 
rate for the automatic weapon 
battery. ELRMA = 20°/sec 


Maximum elevation tracking 
rate for the missile battery. 
ELRMM = 15°/sec 


Maximum Azimuth tracking 
rate for radar controlled gun, 
AZRMG = 12°/sec 


Maximum azimuth tracking 
rate for automatic weapon 
battery. AZRMA = 20°/sec 


Maximum azimuth tracking 
rate for the missile systen. 
AZRMM = 15°/sec 


Maximum effective range for 
radar controlled gun battery. 


Maximum effective range AAW 
battery. RMA = 6,000 feet 


Maximum effective range 
missile system. 
ВММ = 105,600 feet 


Elevation angle for lower 

portion of detection radar 
beam pattern at half power 
point.  THETI = 1.1°2 


Parameter 


ΤΗΕ 


THET 3 


API 


GONT 


IK 


VT 
TRG 


TRM 


ATG 
ΑΤΑ 
ΑΤΗ 
1 


ыр 2 


Ie 


Description 


Elevation angle for detection 
radar, AAW battery. THET 2 =.5 


Elevation angle for detection 
radar, missile battery. 
ТЧЕТ 3 - 1.90 


The constant term in the linear 
equation for slant range versus 
time of slight, gun battery. 

E 


The constant in the expo- 
nential equation for time of 
flight for AAW Battery. 

CONT - ,0002 


Distance between Y grid lines, 
three dimensional terrain. 
YK = 500 feet, 


Target speed, VT = 350 ft/sec 


Crew reaction time, This time 
is measured from time of first 
decision of detection of tar- 
get to track radar lock on, 
This is an average figure. 

TRG = gun system, TRA = AAW 
system, and TRM = missile 
system. TRG 4ο ο ο 

ሬል eC TRM = 30 sec, 


System delay time measured 
from radar lock on to first 
round fired TOF gun battery, 
AAW battery, and missile 
battery. ATG = 20 sec. 
А. ATM - 10 вес. 


Total number of gun batteries. 
TE Ļ 


Total number of AAW batteries 
plusseun batterie. LT 2 = 10 


Total number of gun batteries 
plus AAW batteries plus 
missile batteries. LT 3 = 11 


SIG 


GUNS 


TCRA 


PKSM 


NR 


TIM 


AGUN 


TX 


LY 


YM 


Cross sectional area of target. 
The average value as target 
aspect changes. А = 350 feet 


Angular error of radar: con- 
trolled ‘AA gun when firing at 
2 ከ 1123201710... SIG = 1.5 mille 
Tadians 


The number of guns per gun 
battery. GUNS = 4 


Cyclic rate of fire, AAW 
battery. TCRA = 6 rounds/sec 


Single shot kill probability 
for missile system. PKSM- .^8 


The firing policy of the 
missile battery, NR rounds 
per target. NR= 2 


Time for one round fired from 
the gun battery. TIM = 20 sec 


Number of guns per AAW battery. 
АСОМ = 4 


Memory track period for all 
systems. TX = 15 sec 


‘Index number for the Y coor- 
dinate of the target course 
line terrain points. LY = 1 


Minimum firing range, gun 
battery, RMFG = 500 feet 


Minimum firing range for AA 
battery. RMFA = 500 feet 


Minimum firing range for 
missile system.  RMrM = 9500 
feet ; 


Maximum elevation angle for 
all track radars. FE = 879 


Slope of the linear regression 
equation for time of ء٣٦‎ 
the gun battery. YM = .0007 


2Μ Constant in linear regression 
equation for time of flight 
of missile system. AM = 2,000 


SM Slope of linear regression 
equation. SM = 1200 | 


РС Probability of & kill given ጨ 
hit for the gun battery. This 
is theslithality factors 


PG = «8 

PA Probability of kill given a 
hit for the AAW battery. 
PA = E 


Certain other parameters may be found in the data list. If 
they are not included in the list above, they are not in use 


in the final form of the program and should be disregarded. 


r 


Terrain Data Input. Statement number 7 in the main program 
is the read statement which picks up the data punched on 
84 data cards. Here the variable NTT is set equal to 84, 
Hm ina subsequent use of the program more Terrain points 
were required, changing the value of NTT and providing the 
Same number of data cards provides the change. Statement 
number 8 reads in the coordinates of the Batteries. It is 
essential that the first system type be radar controlled 
gun batteries up to LT 1. Then the next type must be AAW 
batteries up to LT 2 which was defined to be LT 1 Z the 
number of AAW batteries, The last system in order is the 
missile battery. Generally only one missile battery will 
be in range of a single azimuth of approach unless the 


azimuth is in the overlapping zones of fire of two adjacent 
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batteries. Setting the values of LT l, Lt 2, and LT 3 along 
with the proper ordering of the data cards of each type of 
Weapon system constitutes the conversion to a new defensive 


complex. 


Three Dimensional Data Input. If three dimensional data is 
to be used, three new arrays must be introduced for the X, 
Y, and H coordinates of each recorded terrain point. when 
this is done, statement 6 in subroutine SEE must be re- 
placed by à statement using the new altitude array for 


H(J,ly) presently in the program. 
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CHAPTER IV 
PROGRAM VERIFICATION, RESULTS, AND ANALYSIS 


The underlying factor in the model design from concep- 
tion has been that the air defense capability is maximized. 
This means that the attacking aircraft situation is the 
worst possible condition for survival. This fact is appar- 
ent when considering the general model structure presented 
in Chapter II. It should be recalled that the aircraft 
was required to penetrate the defense from out of range to 
the defended area flying a straight and level course, with 
afraid size of one aircraft, and no electronic counter 
measures. In addition, the defense was assumed to be alerted, 
manned, and concentrating on the single intruder. The air- 
craft survivability index calculated in the model when 
submitted to a programing algorithm for minimization within 
the aircraft flight constraints finds the flight path which 
minimized the probability of kill under the worst possible 
attacking conditions. In actual application, anything the 
aircraft may do relieving the sever restrictions, wholly or 
in part, increases the survivability index. The act of 
following the programed output recommended flight path, 
penetrating the defense with more than one aircraft in the 
raid, and using-electronic counter measures will increase 
survivability. These facts must be considered when evalu- 


ating the model's results. 


Ca 


With the development of the model and the program for 
the computer, a verification of its main computing aspects 
is in order, Following this, the computer results will be 


presented and analyzed. 


ጭ 


Verification of Computer Results 


Ihe two main problems solved by the computer within 
the mathmetical model are the calculations for determining 
line of sight, in two dimensional terrain only, and the 
calculations for the probability of kill given a positive 
course line. All other computations within the program 
Support these main requirements. 

Terrain Evaluation. In order to verify the computer solu- 
tion of the line of sight, a sample of computer print out 
c mm in Figures 18 and 19 is presented, Figure 18 is 
a case where the coordinates of the point of maximum radar 
detection, ج‎ 26140. The range of first detection is 
163,060 feet since the coordinates of Battery 6 are given 


as B = 189,200 on the X axis and 6700 feet altitude, The 


- 


LPB ~ 
the terrain point number just greater than Xn 


comparison slope S 00973, EK = 46, and II = 57 "T1 39 


yk? KK is the 


terrain point number just less than В 6° Thus, the program 
has determined the number of intervening terrain points to 


examine for slope comparison. As explained in Chapter II, 


if any slope Sip is less than S PB then line of sight as 


corrected by the distance Z exists. Observing Figure 18 


again, it is seen that 5 562) through S 


L peer) are greater 


L 


ne 


Computer line of Sight Computations 
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Figure 18 
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Figure 19 
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than S When 5 020) is compared to S 15 is found 


LPB' LPB' 

to be less and thus masking line of sight, 2 = 70338.14 
Was calculated from equation (22) Chapter II. A new value 
Dor x 
n 


is given as X = 118862 from equation (23): Using 


K K 
the assumption that targets fly over the terrain when the 
flight path passes through terrain, a new comparison slope 


is calculated and X, is changed from 26140 to 118862, 


k 
Three more terrain obstructions exist such that the final 
value of X yk = 183740, From equation (37), Χρ =-1953.7 
Which indicates that no time is available to fire any 
rounds and battery 6 contributes no probability of kill at 
terrain point 62, and altitude 7. Figure 19 illustrates a 
case where terrain masking occurs once but a distance of 


4830.9 feet for X. remains for firing remembering that the 


D 
time for firing the first and last round have already been 
removed from Хр. Therefore three rounds are fired, one 
more than shown on Figure 19 which is the last round fired. 
The slope comparison calculation is performed 18,490 times 
less those cases where negative elevation angles occur and 
the terrain point is out of range of the battery being 
evaluated, 

Probability Calculations. <A sample of computer output of 
probability calculations is found in Figures 19, 20, and 2l. 
Figure 19 is a sample of a radar controlled gun battery. 

As was seen, three rounds were fired with a probability of 


kill of P (1) = ,46, This is the battery kill probability 


resulting from 12 rounds fired, 3 per gun at a range of 


(2 


Computer Probability Calculations 
Automatic Weapon Battery 
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Figure 20 


Computer Probability Calculations 
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Figure 21 
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5089, 3059, and 1280 feet. Due to the short range, the Р, 
DE and, 2%. Looking st Figure 21, a run of 17 
rounds was fired. Here chesshbiLude eliminates terrain 
masking with a course length of X, = 33762.9 feet. However, 
the range has increased to 23,526 feet so that the Py vary 
Ггоп “00005750 OO The battery kill probability, 

P, Ur) = 6. This probability is the same as for battery 1 
in the earlier example, The short range im tHe other ехалоша 
accounts for the high probability for only 12 rounds fired, 
Figure 20 is the calculations for an automatic weapon battery. 
Fifty rounds were fired in this case. Actually more rounds 
could have been fired but the constraint on this weapon 
system limits firing to 50 rounds. The range is about 4500 


feet on the average with P, varying from .000016 to .0000217. 


1 
The battery kill probability is Б) = „0006, 
Program Results. Appendix V is a set of sample curves wnich 
plot altitude versus cumulative probability of kill for 
specific terrain points. The program prints out the curve 
for each terrain point examined. Included in the wsev are 
lllustrative cS of program results. 

Appendix VI is a set of curves plotting probability 


cost index versus altitude above terrain. This curve is 


obtained by the following relationship; 


= | |- (53) 
P, ur i +Á σοι) 
Ks ο Ν 
IEF (< 30 
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Where 
Poi = the probability that the aircraft pilot 
detects the navigational target. 


Шие values for Ppi were obtained from the computer program 


VISTRAC in reference 1. 


Analysis of Results 

The purpose of the mathmetical and computer models of 
the air defense complex is to compute and plot a curve of 
provability of kill versus altitude at each terrain point. 
Having obtained these curves, it is necessary to analyze 
them in order to determine the reasons for their shape. 
With this information, the application for the model's 
output can then be considered. 

Terrain Effects. Опе of the most significant contributors 
to the shape of these curves is the terrain within which 
the air defense is situated. As was seen in Chapter II, 
Figures 5 through 7, the terrain has different formations 
preducing different curves. The curves in Appendix W 
illustrate this fact. 

Due to inaccessibility, many prominent terrain points 
are not occupied by firing batteries and only cause serious 
masking to batteries located on accessible terrain with 
less altitude than some others. The result is that at low 
altitudes of approach, some batteries may be masked and 
Contribute no probability of kill to the target. whereas 


at higher altitudes the target is unmasked soon enough to 


ЙЭ 


be fired upon by the battery. Figure 22 illustrates this 
concept, The line of sight to altitude 1 intersects tne 
course line at A and Хуу = BE. From this result it 
is clear that the probability of kill at altitude 2, 
since X 


P > ا‎ due to increasing the altitude 


кк2 ? "ga D2 
of approach. With the effect of terrain masking above, 
more rounds will be fired as altitude is increased, 

It can be observed in Appendix V that some curves 
appear to have discrete jumps or regions of discontinuity. 
This is caused by terrain effects as well, These jumps 
occur, when at the lower altitudes, a particular battery 
with high single shot kill probability is masked. At the 
mext altitude of approach considered, the battery is sudden- 
ly unmasked and a large jump in probability is the result. 
Effect of System Parameters. The battery effectiveness 
envelope is another influencing factor to curve shaping. 
This @fivelope was shown in Figure l, Chapter IL, ĪIn.Figūre 
1, notice the curve of the envelope as altitude increases. 
This curve is an assumption which approximates the effect- 
iveness of the average type system in use, If this model 
is used to simulate a weapon system with a radically differ- 
ent envelope, then a modification may be in order. The 
effect of the assumed envelope is also illustrated in 
Figure 22, At altitude 3, the line of sight intersects 
the target course line at G, but is beyond the range of the 
gun system, Therefore, ХЗ = CF since Ç is’ tne point Wire me 


the target course line intersects the weapon system envelope. 
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Figure 22 


It is obvious that CF « BF by the amount of BC' from which 


jt is determined that X = 2 and accordingly P 


рз ^ ^p Kk3 < Ркк: 
This then illustrates why the curves tend to slope down as 
altitude is increased over some regions of altitude. 
Another system parameter wnich effects the output 
curves is radar performance. Ти Chapter II, the half power 
point of a radar beam pattern is assumed to be the con- 
servative limit of usable volume coverage provided,  There- 
fore, any target in this model falling within the beam 


pattern is assumed detected with probability 1. This 


assumes a lateral range curve to be 
ዬ 


ጨ | ይመ "ο 
рса 71-8 (55) 


where 


I^ 


©) KER 
FQ * 1o . в (55) 


-— > m“ 


y 


Substituting (55) into (54) results in 


- l 1 SR 
P(X) = (56) 
These equations are found in reference 11. Figure 23 is a 
plot of such a lateral range curve. The model assumed no 
maximum value for R as a parameter input. The maximum value 
R can take on is limited by the length of the terrain coor- 
dinate system. Maximum R used for tne computer program is 


47 miles. A nominal value for battery acquisition radars is 


65 miles. If a coordinate system were longer than the 


52 


Lateral Range Curve 





Figuro.23 
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maximum rated range of a particular system, then a value of 
В Гог radar detection would nave to be assigned, further, 
ìf the radar detection range normally used in the program 
were at the maximum rated range, the assumption of equation 
(56) may introduce significant error and a different radar 
detection rate would have to be incorporated into the model. 
The reason this factor was not developed was due to the 
extremely low altitudes of target approach where a detection 
range does not become greater than & value much less than 
the maximum rated range of the radar. Such a development 
would constitute a needless complication. 

Two other parameters were used to perform a sensitivity 
analysis to determine the effect of a change in these 
parameters on the probability curves. The minimum and 
maximum range of the AAW battery was originally 500 feet 
and 6000 feet respectively. These values were changed to 
1500 and 9000 feet for one run Of the computer. Curve 
number 8 in Appendix V is the result of using both para- 
meters plotted on one graph? Curve I is the гезити ог tHE 
increased parameters. These two curves are not radically 
different, but the increased range did cause a general 
increase in probability for most values of altitude. The 
rate of fire and total number of rounds fired per run pen 
gun is a more powerful factor in-the model. Curve 5 in 
Appendix V is & curve for AAW batteries firing up to 200 
rounds per gun per run at 1000 rounds per minute with 10 


guns firing independently per battery. The result of this 


a^ 


situation is a very high probability curve which is obvious 
since if sufficient rounds are fired high probabilities will 
result. Practical constraints will not permit this type of 
performance in the usual case. 

Battery ZEuolacement Effects. The probability curves gener- 
ated by the terrain in the computer model are found to be 
date different from terrain point to terrain point, "NEN 
curves from adjacent terrain points are similar in the case 
Where terrain does not change rapidly. Except for this case, 
the curves have different shapes. <A different curve will 
obviously result if a different battery emplacement scheme 
is implemented since new battery to terrain masking situa- 
tions occur. As tne program calculates probabilities for the 
different terrain points, new situations occur such as a 
battery on a high peak is in range or a battery behind a 
peak is unmasked. This variation in terrain conditions is 
why the curves are not exactly related even though the 
probabilities calculated are cumulative, The curves for 
terrain points 47 and 69, curves 3 and 4, Appendix V, when 
compared, are found to be quite different. An additional 
factor to be kept in mind is that the altitudes above 
terrain are not at the same absolute altitude. Altitude 

600 at terrain point 47 is not necessarily the same altitude 
at terrain point 69, There is no intention that there be 
any relationship between the different terrain points. The 


probability curve for each terrain point is separate and 
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distinct from any other unless, of course, ther is little 


terrain change between them except the X coordinate distance. 


Model Application 

In Chapter I, it was pointed out that the work done 
in reference 1 found the relationship between probability of 
radar detection and altitude above terrain to be approxi- 
mately linear. As a result, the computer program incorpor- 
ated a linear fit to the data. Following this the proba- 
bility of pilot navigation was added to the linear function. 
The functional relationship was necessary for use in the 
simplex algorithm for linear program solution for the optimal 
altitude of approach. At the beginning of this work, the 
concept was to develop the air defense model as an input 
to the program in reference l: After examining thevoucpur 
curves of Appendix V a linear fit did not appear feasible 
for the data. In order to determine how bad a linearity 
assumption would be, a regression equation was calculated 
from the data of terrain point 62 in Appendix V. Curve 
number 4 is a plot of the equation. The function derived 


45 given by 


У = .128 + .000369 X (57 


To submit the output data to a linear fit would introđuce 
serious misrepresentation of the model capability. This 
fact is obvious from the results of the linear fit to the 


data of terrain point 62. 
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Probability Cost Index. The curves in Appendix VI are a 
probability cost index. The values are not true proba- 
bilities since they are the sum of two probabilities as 
shown in equation (53). The term cost is defined to be a 
measure of the material effort necessary to penetrate the 
defense, reach the target, and deliver the ordnance. Since 
either one of the events constitutes failure in the mission, 
the higher the value of the cost index the greater the 
chance that the mission will be a failure. Consequently, 
the survivability of the aircraft in the penetration effort 
is maximized when the cost index is minimized. 

When examining the curves in Appendix VI, it is ob- 
served that each curve contains a global minimum. The 
curves then possess a point which minimizes the cost index 
and provides the optimal survivability index number at the 


associated altitude, 
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CHAPTER V 


CONCLUSIONS AND EXTENSIONS 


ThHeVdevewoulent or vie alr defense Model Within the 
assumptions and limitations is complete, As a result of 
the output from the computer model it is concluded: 

(1) That the probability calculations developed a 
result that varies with terrain, air defense systems, and 
employment tactics. Consequently, no assumption of linear- 
ity between probability of kill and altitude above terrain 
is justified. In fact, an investigation of the degree of 
polynomial that a least square fit would require to model 
the™curves in Appendix Vi is necessary; ‘Such Әп ሳሙ 
gation wourd™be required to "discover Wow the"derree would 
vary with terrain points using the same system parameters. 
It would also have to determine how the degree would vary 
with different air defense systems and terrain. Before 
an algorithm can be found to solve for the minimum point 
on the cost curve, the curve itself must be functionally 
0009+ 

(2) That no practical results would be obtained from 
inserting the output from the air defense model Into the 
program in reference 1. If such an attempt were made, tne 
simplex algoritnm used in reference 1 would obtain a solu- 


tion, but it would not be known if the Selutien were the 
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global minimum. Further, the linear regression curve would 
grossly missrepresent the probabilities calculated by the 
air defense model. 

(3) That further investigation 25 necessary to discover 
the feasibility of fitting the curves in Appendix VI to 
some degree polynomial such that the degree chosen is 
acceptable for the curves of all the terrain points. Ifa 
single degree polynomial is not feasible, then a method of 
finding the best fit for each curve over a specified range 
of altitudes must be adopted, 

(4) That having determined the functional relation- 
ship for the probability cost curve, a non ۹99237 
algorithm must be found that solves for the global minimum 
of each curve within the aircraft flight path constraints 


used in reference 1, 


Air Defense Model Extensions 

Two major assumptions in the air defense model should 
be investigated further. The assumption of independence 
between rounds fired successively from the same gun may 
mewroduce serious error. „Везеатецьыисо the determinati onni 
the relative weignts given to correlated and uncorrelated 
error should give an indication of Mow unter isēic Le 
assumption of independence really is. If independence is 
unacceptable then means for finding the conditional proba- 
bilities associated with dependently fired rounds should 


be studied, 
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An additional problem for further study is the assump- 
tion that the gun fire snot pattern is circular normal. 

Much evidence points to the fact that the shot pattern is 
actually elliptical normal. Empirical equations for the 
elliptical normal error pattern such as equation (42) and 
(48), Chapter II in the circular normal case should be used 
if such equations exist, 

The final extension is in the area of three dimensional 
terrain. As Was previously БЕПОЯ out, the Fortran code 
written for evaluation of line of sight in three dimensional 
terrain was not verified. With the availability of digital- 
ized terrain this portion of the program should be verified. 
However, one serious disadvantage in the concept used in 
this model is the requirement that the target course Lime 
be parallel to the X axis of the coordinate system. The 
problem lies in the necessity of assigning new altitude 
values for each different route of approach into the target. 
Another method for determining line of sight between 2 
points in 3 dimensional terrain is found in reference 6, 
This method requires a least square fit to the terrain 
points lying on the line between the two points for wnich 
line oP Sight 3e being determined, Then the method deter- 


mines if the curve extends above the line of sight. 
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TERRAIN AND BATTERY COORDINATES 





COORDINATES OF MAJOR TERRAIN POINTS 
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APPENDIX III 


GLOSSARY OF FORTRAN VARIABLE NAMES 


The following is a list of Fortran variable names in 


the order of appearance. This list excludes the system 


parameters which were presented in Chapter III. 


КК, в) 


H(K,LY) 


SIKE) 


BX(J) 


BY(J) 


BH(J) 


RANGE(I) 


АН(Т) 


DELTX 


DELTY 


DELTH 


THE X COORDINATE OF THE KTH TERRAIN 

POINT WITH THE LY TH Y COORDINATE INDEX 
THE H COORDINATE OR ALTITUDE OF THE KTH 
TERRAIN POINT WITH LY TH Y COORDINATE 
INDEX 

THiS COORDINATE OF THHeHTH TERRAIN POT 
WITH LY TH Y COORDINATE INDEX, SET EQUAL 
ΤΟ O FOR TWO DIMENSIONAL TERRAIN 

X COORDINATE OF THE JTH BATTERY 

Y COORDINATE OF THE JTH BATTERY SET EQUAL 
TO O FOR 2 DIMENSIONAL TERRAIN 

H COORDINATE OF THE JTH BATTERY 

AN ARRAY USED WITH SUBROUTINE UTPLOT 
ALTITUDE ABOVE TERRAIN FOR THE ITH POINT 
THE DIFFERENCE IN X COORDINATES OF THE 
KTH TERRAIN POINT AND THE JTH BATTERY 
THE DIFFERENCE IN Y COORDINATES OF THE 
KTH TERRAIN POINT AND JTH BATTERY 

THE DIFFERENCE IN THE H COORDINATES OF 


THE KTH TERRAIN POINT AND THE JTH BATTERY 
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ELRX 


AZRX 


RM 


THETA 


АТ 


ER 


RS 


፲፻ 


RXS 


ШЕТ 


RSA 


RSE 


MAXIMUM ELEVATION TRACKING RATE FOR THE 
SPECIFIC BATTERY BEING EXAMINED 

MAXIMUM AZIMUTH TRACKING RATE OF THE 
SPECIFIC BATTERY BEING EXAMINED 

MINIMUM RANGE OF THE SPECIFIC BATTERY 
BEING EXAMINED 

MAXIMUM RANGE OF THE BATTERY BEING 
EXAMINED 

DETECTION RADAR LOWER EDGE OF BEAM 
PATTERN ELEVATION ANGL& IN DEGREES OF 
THE BATTERY BEING EXAMINED 

SYSTEM DELAY TIME FOR THE BATTERY BEING 
EXAMINED 

CREW REACTION TIME FOR THE BATTERY BEING 
EXAMINED 

SLANT RANGE TO THE (K,I)TH POINT 
HORIZONTAL COMPONENT OF MINIMUM RANGE 
HORIZONTAL COMPONENT OF RS 

HORIZONTAL COMPONENT OF THE RANGE OP THE 
POINT WHERE THE LINE DRAANN FROM THE 
MAXIMUM ELEVATION ANGLE INTERSECTS THE 
TARGET COURSE LINE 

THE SLANT RANGE WHERE MAXIMUM AZIMUTH 
RATE OCCURS 

THE SLANT RANGE WHERE MAXIMUM ELEVATION 


RATE OCCURS 


2. 


XD 


ROF 


NS ,XNR 


RH 


ТЕ 
GUNS, NGUNS 


R(II) 
SIG(II) 


SIGM 


EP(II) 


TA 


PD(J) 


τα 


NT, XNR 


LE 


THE DISTANCE ON THE TARGET COURSE LINE 
AVAILABLE FOR FIRING ALL BUT FIRST AND 
LAST ROUNDS 

OPEN FORE SLANT RANGE 

TIME AVAILABLE FOR FIRING EXCLUDING 
FIRST AND LAST ROUNDS 

NUMBER OF ROUNDS FIRED PER RUN 
HORIZONTAL RANGE TO MAXIMUM RADAR 
DETECTION POINT ON COURSE LINE 

TIME OF FLIGHT 

NUMBER OF GUNS PER BATTERY, RADAR CON- 
TROLLED GUN BATTERY 

RANGE OF II TH ROUND FIRED 

THE FIRING ERROR RADIUS OF THE II TH 
ROUND FIRED 

FIRING ERROR OF THE RADAR CONTROLLED GUN 
THE PROBABILITY OF KILL OF THE II TH 
ROUND FIRED 

NS 

TT (1-EP(II)) FOR THE RADAR CONTROLLED 
L=I 

BATTERY 

PROBABILITY OF KILL OF THE JTH BATTERY 
NS 

TT (1-EP(II)) FOR AAW BATTERY 

L=I 

NUMBER OF MISSILES FIRED 

NFU 

TT (1-PK(J)) THE PROBABILITY THAT NO 


I= I 
ROUND KILLS THE TARGET 
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PHI(K,I) 


RAM 


ANK 


Rd 


SLPD 


SLP(M) 


ADR 


ET 


PROBABILITY THE AIRCRAFT IS KILLED AT THE 
(K,I)TH POINT IN SPACE 

SUBROUTINE LOS 
HORIZONTAL COMPONENT OF MAXIMUM FIRING 
RANGE 
X COORDINATE ON TARGET COURSE LINE OF 
FIRST RADAR DETZCTION POINT 
X COORDINATE INDEX OF LAST TERRAIN POINT 
BETWEEN POINT IN SPACE FOR CALCULATION 
AND THE BATTERY 
X COORDINATE INDEX OF FIRST TERRAIN POINT 
BETWEEN POINT IN SPACE FOR CALCULATION 
AND THE BATTERY 
THE SLOPE OF THE LINE BETWEEN THE POINT 
OF FIRST RADAR DETECTION AND THE BATTERY 
ΤΗΕ SLOPE OF EACH OF THE M TERRAIN POINTS 
THAT MAY CAUSE TERRAIN MASKING ON EACH 
FIRING RUN 
DISTANCE ON TARGET COURSE LINE FROM POINT 
OF TERRAIN UNMASKING TO THE BATTERY 
THE DISTANCE BETWEEN POINT OF FIRST 
RADAR DETECTION AND OPEN FIRE POINT, ALSO 
CALLED WARNING TIME DISTANCE 
TOTAL DELAY DISTANCE ON TARGET COURSE 


LINE SUBROUTINE SEE 


22 


INDICATION OF WHETHER LINE OF SIGHT EXISTS 
FOR EACH CALLING OF THE SUBROUTINE SEE 
EYE 0 INDICATES LINE OF SIGHT EXISTS 
EYE = 1.0 INDICATES LINE OF SIGHT DOES 


NOT EXIST 
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APPENDIX IV 


FORTRAN PROGRAM 





17751755 


DATE x 623328 


MAIN 


1, MOD 2 


LEVEL 
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APPENDIX VI 


COST VERSUS ALTITUDE FOR SELECTED 
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DATA POINTS FOR THE CURVES IN APPENDICES 
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